Generalized model for dynamic percolation 
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We study the dynamics of a carrier, which performs a biased motion under the influence of 
an external field E, in an environment which is modeled by dynamic percolation and created 
by hard-core particles. The particles move randomly on a simple cubic lattice, constrained by 
hard-core exclusion, and they spontaneously annihilate and re-appear at some prescribed rates. 
Using decoupling of the third-order correlation functions into the product of the pairwise carrier- 
particle correlations we determine the density profiles of the "environment" particles, as seen 
from the stationary moving carrier, and calculate its terminal velocity, V c , as the function of the 
applied field and other system parameters. We find that for sufficiently small driving forces the 
force exerted on the carrier by the "environment" particles shows a viscous-like behavior. An 
analog Stokes formula for such dynamic percolative environments and the corresponding friction 
coefficient are derived. We show that the density profile of the environment particles is strongly 
inhomogeneous: In front of the stationary moving carrier the density is higher than the average 
density, p s , and approaches the average value as an exponential function of the distance from 
the carrier. Past the carrier the local density is lower than p s and the relaxation towards p s may 
proceed differently depending on whether the particles number is or is not explicitly conserved. 

PACS numbers: 05.40. +j, 05.60,+w, 02.50.+S, 05.70.Ln, 47.40. Nm 



I. INTRODUCTION. 

The percolation concept has turned out very useful for understanding transport and 
conduction processes in a wide range of disordered media, as exemplified by ionic conduction 
in polymeric, amorphous or glassy ceramic electrolytes, diffusion in biological tissues and 
permeability of disordered membranes [|]-|3|]. 

Most of the situations discussed in Refs. |l|-|3| pertain, however, to systems with "frozen" 
disorder; that is, the random environment in which a given transport process takes place does 
not change in time. This is certainly the case in many instances, but it is not true in general. 
As a matter of fact, there are many experimental systems in which the static percolation 
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picture does not apply since the structure of the host material undergoes essential structural 
reorganizations on a time scale comparable to that at which the transport itself occurs. A few 
stray examples of such systems include certain biomembranes ||], solid protonic conductors 
||, oil-continuous microemulsions and polymer electrolytes flOHT2 . 



More specifically, ionic transport across a biomembrane, such as, e.g. gramicidin- A, oc- 
curs by the motion of ions through molecular channels along which they encounter potential 
barriers that fluctuate in time. The fluctuations of potential barriers may hinder signifi- 
cantly the transport and constitute an important transport-controlling factor ||. In the 
case of protonic conduction by the Grotthus mechanism || , site-to-site motion of carriers 
occurs only between those neighboring H 2 or NH 3 groups that have a favorable relative 
orientation; thermally activated rotation of these groups is the structural host-reorganization 
process interacting with the carrier motion. Similarly, within oil-continuous microemulsions, 
the charge transport proceeds by charge being transfered from one water globule to another, 
as globules approach each other in their Brownian motion 0-0. Lastly, in polymer elec- 
trolytes, such as, e.g. polyethylene oxide complexed nonstoichiometrically with the ionic 
salt NaSCN, the Na + ions are largely tetrahedrally coordinated by polyether oxygenes, 
but at the same time that Na + ions hop from one fourfold coordination site to another, the 
oxygens themselves, along with the polymeric backbone, undergo large-amplitude wagging 



and even diffusive motion [ 10 12 



Clearly, all the above mentioned examples involve two characteristic time scales, one 
which describes the typical time r between two successive hops of the carrier, while the 
other is associated with a typical renewal time r* of the environment itself; namely, the 
time needed for the host medium to re-organize itself and thereby provide a new set of 
available pathways for transport. Consequently, the static percolation picture applies only 
when the characteristic time r* gets infinitely large. For a finite r* dynamic percolation 
has to be considered, and one encounters quite a different behavior when compared to the 
random environments with quenched disorder. As a result, one observes Ohmic-type or 
Stokes-type linear velocity-force relation for the carrier's terminal velocities as a function of 



the applied field, in contrast to the threshold behavior predicted by the static percolation 
theory. The prefactor in the linear velocity-force relation may depend, however, in a non- 
trivial way on the system's parameters and this dependence consitutes the main challenge 
for the theoretical analysis here. On the other hand, we note that in the above mentioned 
examples of the dynamic percolative environments quite different physical processes are 
responsible for the time evolution of the host medium. Consequently, one expects that the 
prefactor in the Stokes-type velocity-force relation should also be dependent on the precise 
mechanism which underlies the temporal re-organization of the environment. 

Theoretical modelling of dynamic percolative environments has followed several avenues, 
which differ mostly in how the time evolution of the disorder is constrained; Namely, is it con- 
strained (a) by conservation laws or (b) by spatial and temporal correlations in the renewal 
events? Early models of dynamic percolation |~3],[14] described the random environment 
within the framework of a standard bond-percolation model, in which the strength of each 
bond fluctuates in time between zero and a finite value. The dynamics of the host medium 
in these models [ |T3| , [H| was accounted for by a series of instantaneous renewal events. These 



events were assumed to occur at random times, chosen from a renewal time distribution. 
In the renewal process the positions of all unblocked bonds are being reassigned, such that 
after each renewal event a carrier sees a newly defined network. This approach is thus char- 
acterized by a global dynamical disorder without global conservation laws and correlations, 
since the entire set of random hopping rates is simultaneously renewed independently of 
the previous history. Another model characterized by a local dynamical disorder has been 
proposed in Refs. |15| and [fTE], and subsequently generalized to the non-Markovian case 



in Ref. [17]]. This model appears to be similar to the previous one, except that here the 



hopping rates at different sites fluctuate independently of each other. That is, individual 
bonds, rather than the whole lattice change in the renewal events. To describe the dynami- 
cal behavior in the local dynamical disorder case, a dynamical mean-field theory has been 
proposed P3|,T6], based on the effective medium approximation introduced for the analysis 



of random walks on lattices with static disorder ||18|| , and has been generalized to include 



the possibility of multistate transformations of the dynamically random medium [fL9| . More 
recently, several exactly solvable one-dimensional models with global and local dynamical 
disorder have been discussed P0[j . 

In the second approach, which emerged within the context of the ionic conductivity 
in superionic solids, the dynamical percolative environment has been considered as a mul- 
ticomponent mixture of mobile species in which one or several neutral components block 
the carrier component ||21|| . In particular, such a situation can be observed in a superionic 



conductor /3"-alumina, doped with two different ionic species (e.g. Na + and Ba 2+ ), where 
small Na + ions are rather mobile, while the larger Ba 2+ ions move essentially slower and 
temporarily block the Na + ions. Contrary to the previous line of thought, the dynamics of 
such a percolative environment has essential correlations, generated by hard-core exclusion 
interactions between the species involved, and moreover, it obeys the conservation law - the 
total number of the particles involved is conserved. In Ref. [^TJ, the frequency-dependent 



ionic conductivity of the light species has been analysed combining a continuous time ran- 
dom walk approach for the dynamical problem with an effective medium approximation 
describing the frozen environment of slow species. Next, as an explanation of the sharp 
increase of electrical conductivity transition in water-in-oil microemulsions when the vol- 
ume fraction of water is increased towards a certain threshold value, in Refs. @ and f| 
it has been proposed that the charge carriers are not trapped in the finite water clusters, 
but rather a charge on a water globule can propagate by either hopping to a neighboring 
globule, when they approach each other, or via the diffusion of the host globule itself. This 



picture has been interpreted in terms of a model similar to that employed in Ref. [21], with 



the only difference being that here the "blockers" of Ref. 0] play the role of the transient 



charge carriers. In the model of Refs. JF| and ||, in which the host dynamics is influenced 
by spatial correlations and conservation of the number of the water globules involved, the 
conductivity depends hence, on the rate of cluster rearrangement. Lastly, a similar problem 
of a carrier diffusion in an environment created by mobile hard-core lattice-gas particles has 
been analysed in Ref. |22| by using the developed dynamic bond percolation theory of Refs. 



g and H . 



In this paper we propose a generalized model of dynamic percolation which shares com- 
mon features with both bond-fluctuating models of Refs. |13|- |iT| , p||20|| as well as models 



involving mobile blockers of Refs. |2l|]2!q| . The system we consider consists of a host lattice, 
which here is a regular cubic lattice whose sites support at most a single occupancy, hard- 
core "environment" particles, and a single hard-core carrier particle. The "environment" 
particles move on the lattice by performing a random hopping between the neighboring 
lattice sites, which is constrained by the hard-core interactions, and may disappear from 
and re-appear (renewal processes) on the empty sites of the lattice with some prescribed 
ratesQ. In turn, the carrier particle is always present on the lattice, i.e., it can not disappear 
spontaneously, and is subject to a constant external force E. Hence, the carrier performs 
a biased random walk, which is constrained by the hard-core interactions with the "envi- 
ronment" particles, and probes the response of the percolative environment to the internal 
perturbancy or, in other words, the frictional properties of such a dynamical environment. 

An important aspect of our model, which makes it different to the previously proposed 
models of dynamic percolation, is that we include the hard-core interaction between "envi- 
ronment" particles and the carrier molecule, such that the latter may influence the dynamics 
of the environment. This results, as we proceed to show, in the emergence of complicated 
density profiles of the "environment" particles around the carrier. These profiles, as well 
as the terminal velocity V c of the carrier, are determined here explicitly, in terms of an 



*We hasten to remark that diffusive processes, of course, also result in a certain renewal of the 
environment; diffusive processes, as compared to the spontaneous creation and annihilation of 
particles, have however completelely different underlying physics and influence in a completely 
different fashion the evolution of the system, as we proceed to show. Following the terminology of 



Refs. [13 -17, 19, 20 1, we thus choose here to distinguish between diffusive and creation/annihilation 



processes, referring to the latter as the renewal ones. 
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approximate approach of Ref . p3 , which is based on the decoupling of the carrier-particle- 



particle correlation functions into the product of pair-wise correlations. We show that the 
"environment" particles tend to accumulate in front of the driven carrier creating a sort of 
a "traffic jam", which impedes its motion. Thus the density profiles around the carrier are 
highly asymmetric: the local density of the "environment" particles in front of the carrier 
is higher than the average and approaches the average value as an exponential function of 
the distance from the carrier. The characteristic length and the amplitude of the density 
relaxation function are calculated explicitly. On the other hand, past the carrier the local 
density is lower than the average: We show that depending on the condition whether the 
number of particles in the percolative environment is explicitly conserved or not, the local 
density past the carrier may tend to the average value either as an exponential or even as 
an algebraic function of the distance, revealing in the latter case especially strong memory 
effects and strong correlations between the particle distribution in the environment and the 
carrier position. Further on, we find that the terminal velocity of the carrier particle de- 
pends explicitly on the excess density in the "jammed" region in front of the carrier, as well 
as on the "environment" particles density past the carrier. Both, in turn, are dependent 
on the magnitude of the velocity, as well as on the rate of the renewal processes and the 
rate at which the "environment" particles can diffuse away from the carrier. The interplay 
between the jamming effect of the environment, produced by the carrier particle, and the 
rate of its homogenization due to diffusive smoothening and renewal processes, manifests 
itself as a medium-induced frictional force exerted on the carrier, whose magnitude depends 
on the carrier velocity. As a consequence of such a non-linear coupling, in the general case, 
(i.e. for arbitrary rates of the renewal and diffusive processes), V c can be found only im- 
plicitly, as the solution of a non-linear equation relating V c to the system parameters. This 
equation simplifies considerably in the limit of small applied external fields E and we find 
that the force-velocity relation to the field becomes linear. This implies that the frictional 
force exerted on the carrier particle by the environment is viscous. This linear force- velocity 
relation can be therefore interpreted as the analog of the Stokes formula for the dynamic 
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percolative environment under study; in this case, the carrier velocity is calculated explic- 
itly as well as the corresponding friction coefficient. In turn, this enable us to estimate 
the self-diffusion coefficient of the carrier in absence of external field; we show that when 
only diffusive re-arrangement of the percolative environment is allowed, while the renewal 
processes are suppressed, the general expression for the diffusion coefficient reduces to the 
one obtained previously in Refs. |24|] and p5 |. We note that the result of Refs. |24f and [25 



is known to serve as a very good approximation for the self- diffusion coefficient in hard-core 
lattice-gases |^6| . 



We finally remark, that a qualitatively similar physical effect was predicted recently for 
a different model system involving a charged particle moving at a constant speed a small 
distance above the surface of an incompressible, infinitely deep liquid. It has been shown 



in Refs. p7| , |28[| , that the interactions between the moving particle and the fluid molecules 
induce an effective frictional force exerted on the particle, producing a local distortion of the 
liquid interface, - a bump, which travels together with the particle and increases effectively 
its mass. The mass of the bump, which is analogous to the jammed region appearing in our 
model, depends itself on the particle's velocity resulting in a non-linear coupling between 
the medium-induced frictional force exerted on the particle and its velocity |27],|23. 

The paper is structured as follows: In Section II we formulate the model and introduce 
basic notations. In Section III we write down the dynamical equations which govern the 
time evolution of the "environment" particles and of the carrier. Section IV is devoted to 
the analytical solution of these evolution equations in the limit t — > oo ; here we also present 
some general results on the shape of the density profiles around stationary moving carrier and 
on the carrier terminal velocity, which is given implicitly, as the solution of a transcendental 
equation defining the general force- velocity relation for the dynamic percolative environment 
under study. In Section V we derive explicit asymptotic results for the carrier terminal 
velocity in the limit of small applied external fields E and obtain the analog of the Stokes 
formula for such a percolative environment; as well, we present here explicit results for the 
friction coefficient of the host medium and for the self- diffusion coefficient of the carrier in the 



absence of external field. Asymptotic behavior of the density profiles of the "environment" 
particles around the carrier is discussed in Section VI. Finally, we conclude in Section VII 
with a brief summary and discussion of our results. 



II. THE MODEL. 

The model for dynamic percolation we study here consists of a three-dimensional simple 
cubic lattice of spacing a, the sites of which are partially occupied by identical hard-core 
"environment" particles and a single, hard-core, carrier particle (see Fig.l). For both types 
of particles the hard-core interactions prevent multiple occupancy of the lattice sites; that is, 
no two "environment" particles or the "carrier" and an "environment" particle can occupy 
simultaneously the same site, and particles can not pass through each other. 

The occupation of the lattice sites by the "environment" particles is characterized by the 
time-dependent occupation variable r/(r), r being the lattice- vector of the site in question. 
This variable assumes two values: 

!1, if the site f is occupied 
(1) 
0, if the site r is empty 

Next, we assume the following dynamics of the "environment" particles: The particles can 
spontaneously disappear from the lattice, and may re-appear at random positions and ran- 
dom time moments, which is reminiscent of the host medium dynamics stipulated in Refs. 
13| 17, |19 , p0 . We refer to these two processes generally as renewal processes. In addition, 



the environment particles move randomly within the lattice by performing nearest-neighbor 
random walks constrained by the hard-core interactions, which is the main feature of the 
approach in Refs. |f2T| , |22|| . We stipulate that any of the "environment" particles waits a time 
5t, which has an exponential probability distribution with a mean r*, and then chooses from 
a few possibilities: (a) disappearing from the lattice at rate g, which is realized instanta- 
neously, or (b) attempting to hop, at rate 1/6, onto one of 6 neighboring sites. The hop is 
actually fulfilled if the target site is not occupied at this time moment by any other particle; 
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otherwise, the particle attempting to hop remains at its initial position, and (c) particles 
may re-appear on any vacant lattice site with rate /. 

Note that, for simplicity, we assumed that the characteristic diffusion time and the 
renewal times of the "environment" particles are equal to each other. These times, i.e. 
Tdif, mean creation time r cr and mean annihilation time r an may, however, be different, and 
can be restored in our final results by a mere replacement I — > It* /r^f, f — > fr*/r cr and 

g -> gr*/r an . 

Note also that the number of particles is not explicitly conserved in such a dynamical 
model of the environment, which happens because of the presence of the renewal processes; 
the particles diffusion, on contrary, conserves the particles number. However, in the absence 
of attractive particle-particle interactions and external perturbances, the particles distribu- 



tion on the lattice is uniform and the average occupation p(t) = 77(f) of the lattice tends, 
as t — > oo, to a constant value, p s = //(/ + g). This relation can be thought of as the 



Langmuir adsorption isotherm [55] 



Hence, the limit r^/ — > oo (or, / — > 0) corresponds to the ordinary site percolation model 
with immobile blocked sites. The limit f,g—>0, (r cr , r an — > oo), while keeping the ratio f/g 
fixed, f/g = p s /(l — p s ), corresponds to the usual hard-core lattice-gas with the number of 
particles conserved. 

At time t = we introduce at the origin of the lattice an extra particle, the carrier, 
whose role is to probe the response of the environment modeled by dynamic percolation to 
an external perturbance. We stipulate that only the carrier out of all participating particles 
can not disappear from the system, and moreover, its motion is biased by some external 
constant force. As a physical realization, we envisage that the carrier is charged, while all 
other particles are neutral, and the system is exposed to constant external electric field E. 
The dynamics of the carrier particle is defined as follows: We suppose that the waiting 
time between successive jumps of the carrier has also an exponential distribution with a 
mean value r, which may in general be different from the corresponding waiting time of 
the environment particles. Attempting to hop, the carrier first chooses a hop direction with 
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probabilities 



p^ = exp 



^•e,)]/E^[f(£-4)], (2) 



2 



2 



where (3 is the reciprocal temperature, e u (or e M ) stand for six unit lattice vectors, v, fi = 
{±1, ±2, ±3}, connecting the carrier position with 6 neighboring lattice sites, and (E ■ e v ) 
denotes the scalar product. We adopt the convention that ±1 corresponds to ±X, ±2 
corresponds to ±Y while ±3 stands for ±Z. The jump is actually fulfilled when the target 
lattice site is vacant. Otherwise, as mentioned the carrier remains at its position. For 
simplicity we assume in what follows that the external field is oriented along the X-axis in 
the positive direction, such that E = (E, 0, 0). Note also that for the choice of the transition 
probabilities as in Eq.(Q), the detailed balance is naturally preserved. 

III. EVOLUTION EQUATIONS 

Let P(R c ,rj;t) denote the joint probability that at time moment t the carrier occupies 
position R c and all "environment" particles are in configuration i] = {f](r}}- Next, let rf' 11 
denote particles' configuration obtained from rj by exchanging the occupation variables of 
the sites r and r + ep, i.e. r](f) <-> rj(f+ ep), and ff be the configuration obtained from t] by 
changing the occupation of the site f as rj(r) <-> 1 — r](r). Clearly, the first type of process 
appears due to random hops of the "environment" particles, while the second one stems from 
the renewal processes, i.e. random creation and annihilation of the "environment" particles. 
Then, summing up all possible events which can result in the configuration (R c , rf) or change 
this configuration for any other, we find that the temporal evolution of the system under 
study is governed by the following master equation: 

d t P(R c , m t) = -±-Yl E { p (^, rf^- 1) - P(R C , V ;t)} + 

A* r^R c -e p ,R c 

+ ~ E?V {(! - v(Rc)) P(Rc - e n , rj; t) - (l - V (R C + P(R CJ m t)} + 
T v 



+ ^ E {(i-v(^)P(Rc,tf;t)- v (r>)P(R c , m t)} 

r^Rc 
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+ {v(^)P(Rcf;t)-(i- v (^)P(R c , m t)}. (3) 

Note that the terms in the first (resp. second) line of Eq.(||) describe random hopping 
motion of the "environment" particles (resp. biased motion of the carrier) in terms of the 
Kawasaki-type particle-vacancy exchanges, while the terms in the third and the fourth lines 
account for the Glauber-type decay and creation of the "environment" particles. 

A. Mean velocity of the carrier and correlation functions. 

From Eq.(H) we can readily compute the velocity of the carrier. Multiplying both sides 
of Eq.(Q) by (R c ■ ei) and summing over all possible configurations (R c ,v) we find that the 
carrier's mean velocity V c (t), defined as 

V e (t) = i(5^), (4) 

obeys: 

V c (t) = - {pj(l - fc(e x ; t)) - p_ x (l - fe(e_ i; t))} , (5) 

T 

where k(\;t) stands for the carrier-" environment" particles pair correlation function 

k{X;t) = J2r ] (R c + \)P(R c ,r l ;t). (6) 

Rc,v 

In other words, k(X;t) can be thought of as the density distribution of the "environment" 
particles as seen from the carrier which moves with velocity V c (t). 

Hence, V c (t) depends explicitly on the local density of the "environment" particles in the 
immediate vicinity of the carrier. Note that if the "environment" is perfectly homogeneous, 
i.e. if for any A the density profile is constant, k(X;t) = p s , which immediately implies 
decoupling between rj(R c + A) and P(R C , V> t) m Eq.(ff), then we obtain from Eq.(||) a trivial 
mean-field-type result 

K (0) = (Pi-P-i)(l-P*A (7) 

T 
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which states that the frequency of jumps of the carrier particles (t _1 ) only gets renormalized 
by a factor 1 — p s , which gives the fraction of successful jumps. 

The salient feature of our model is that there are essential backflow effects. The carrier 
effectively perturbs the spatial distribution of the "environment" particles so that stationary 
density profiles emerge. This can be contrasted to the earlier dynamic percolation models 
|l3|- |l7i|l9H 22|] in which the carrier had no impact on the embedding medium and hence there 
was no re-arrengement of the host medium around the carrier particle. As a consequence 
k(X; t) ^ p s , and k(X; t) approaches p s only at infinite separations from the carrier, i.e. when 
|A| — ► oo. Therefore, we rewrite Eq.(||) in the form 



which shows explicitly the deviation of the mean velocity of the carrier from the mean-field- 



From Eq.(|]) it follows that in order to obtain V c (t), it suffices to compute k(e±i,t). 
Consequently, we have to evaluate the equation governing the time evolution of the pair 



V c (t) = K (0) - -{pi{k{e l -t)-p s )- V -i{p s - fc(e_ i; t))}, 



(8) 



type result in Eq.(|7]) due to the formation of the density profiles. 



B. Evolution equations of the pair correlation functions. 



correlation functions. Multiplying both sides of Eq.(|3|) by rj(R c ) and summing over all 
configurations (R c ,rj), we find that k(X;t) obeys 



WA; t) = ± £(V„ - 5 lg V^)k(X; t) - { -^-k(X; t) + I + 




where V M denotes the ascending finite difference operator of the form 



v M /(X) = /(X + ^)-/(X), 



(10) 



and 




0, otherwise. 



1, if the site f— f 



(11) 
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The Kroneker-delta term 5% g signifies that the evolution of the pair correlations, Eq.(^), 
proceeds differently at large separations and at the immediate vicinity of the carrier. This 
stems from the asymmetric hopping rules of the carrier particle defined by Eq.(|2|). 

Note next that the contribution in the second line in Eq.(|]), which is associated with 
the biased diffusion of the carrier, appears to be non-linear with respect to the occupation 
numbers, such that the pair correlation function gets effectively coupled to the evolution of 
the third-order correlations of the form 

T(A, e u ; t) = £ rj{R c + X) V (R C + e u )P(R c , 77; t). (12) 

Rc,v 

That is, Eq. (|D is not closed with respect to the pair correlations but rather represents a first 
equation in the infinite hierachy of coupled equations for higher-order correlation functions. 
One faces, therefore, the problem of solving an infinite hierarchy of coupled differential 
equations and needs to resort to an approximate closure scheme. 



C. Decoupling Approximation 

Here we resort to the simplest non-trivial closure approximation, based on the decoupling 
of the third-order correlation functions into the product of pair correlations. More precisely, 
we assume that for A 7^ e„, the third-order correlation fulfils 

R C ,T] 

J2v(Rc + \)P(RciV,t)\ l^2rj(R c + e v )P(R c ,rj;t)Y (13) 
\Rc,v > \Rc,v / 

or, in other words, 

J2 V(R C + A)t?(4 + e u )P(R c , 77; t) « fc(A; t)k{e v] t) (14) 
Rc,v 

The approximate closure in Eq. (0) has been already employed for studying related models 
of biased carrier diffusion in hard-core lattice gases and has been shown to provide quite an 
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accurate description of both the dynamical and stationary-state behavior. The decoupling 
in Eq.(|i~4|) was first introduced in Ref. J^J to determine the properties of a driven carrier 
diffusion in a one-dimensional hard-core lattice gas with a conserved number of particles, 
i.e. without an exchange of particles with the reservoir. Extensive numerical simulations 
performed in Ref. |23| have demonstrated that such a decoupling is quite a plausible ap- 
proximation for the model under study. Moreover, rigorous probabilistic analysis of Ref. 
|30| has shown that for this model the results based on the decoupling scheme in Eq.(|l4|) 



are exact. Furthermore, the same closure procedure has been recently applied to study 
spreading of a hard-core lattice gas from a reservoir attached to one of the lattice sites [|3T] . 



Again, a very good agreement between the analytical results and the numerical data has 
been found. Next, the decoupling in Eq.([T4]) has been used in a recent analysis of a biased 
carrier dynamics in a one-dimensional model of an adsorbed monolayer in contact with a 
vapour phase [Q, i.e. a one- dimensional version of the model to be studied here. Also in 
this case an excellent agreement has been observed between the analytical predictions and 
the Monte Carlo simulations data |32[| . We now show that the approximate closure of the hi- 
erarchy of the evolution equations in Eq. (|TJ]) allows us to reproduce in the limit /, g = and 
/ / g = const the results of Refs. [23]| and which are known (see e.g. Ref. [|2B|]) to provide 
a very good approximation for the carrier diffusion coefficient in three-dimensional hard-core 
lattice gases with arbitrary particle density. We expect therefore that such a closure scheme 
will render a plausible description of the carrier dynamics in a three-dimensional generalized 
dynamic percolation model. We base our further analysis on this approximation. 

Making use of Eq. fll4]) , we find from Eq.(|9]) that the pair correlations obey the following 
equations: 

d t k(\;t) = -^Lk(\-,t)+^ (15) 
or* t* 

which hold for all lattice sites except for those at the immediate vicinity of the carrier, i.e. 
for all A except for A = {0, e±\, e± 2 , e± 3 }, while at the sites adjacent to the carrier one has 

d M^ t) = ^ (L + A v {t))k{e v - 1) + (16) 
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where v = {±1, ±2, ±3}. The operators L and coefficients A u (t) are given explicitly by 

L = $>,(t) V, - (17) 

and 

6r* 

A„(t) = 1 + — - A;(e M ; *)), (18) 
tr 

where V M has been defined previously in Eq. (|TI1|) , /z = {±1,±2,±3}. It is important to 
emphasize that all coefficients A^t) = A^(E,V c ;t), i.e. are functions of both the applied 
field and the carrier velocity. 

Now, several comments about equations (|i~5|) and ( |16|) are in order. First of all, let us 
note that Eq. (|i~6|) represents, from the mathematical point of view, the boundary conditions 
for the general evolution equation fllB] ), imposed on the sites in the immediate vicinity of 
the carrier. Equations (|I5|) and ([16]) have a different form since in the immediate vicinity 
of the carrier its asymmetric hopping rules perturb essentially the "environment" particles 



dynamics. Equations (|15D and (|16|) possess some intrinsic symmetries and hence the number 
of independent parameters can be reduced. Namely, reversing the field, i.e. changing E — > 
—E, leads to the mere replacement of k{e\]t) by /c(e_i;t) but does not affect k(e u ;t) with 
v = {±2, ±3}, which implies that 

h(e 1 ;t)(-E) = k(e- 1 ;t)(E), and k(e„; t)(-E) = k(e v ; t)(E) for z/ = {±2,±3}, (19) 

Besides, since the transition probabilities in Eq.(Q) obey 

P2 = P-2 =P3= P-3 (20) 

one evidently has that 

k(e 2 ; t) = fc(e„ 2 ; t) = k(e 3 ; t) = fc(e_ 3 ; t), (21) 

and, by symmetry, 

A 2 (t) = A_ 2 (t) = A 3 (t) = A- 3 (t) (22) 
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which somewhat simplifies equations (|T5|) and flUf). Lastly, we note that despite the fact 
that using the decoupling scheme in Eq.(Jl^) we effectively close the system of equations 
on the level of the pair correlations, the solution of Eqs. (|15D and ([16]) still poses serious 
technical difficulties. Namely, these equations are strongly non-linear with respect to the 
carrier velocity, which introduces the gradient term on the rhs of the evolution equations 
for the pair correlation, and depends by itself on the values of the "environment" particles 
densities in the immediate vicinity of the carrier. Below we discuss a solution to this non- 
linear problem, focusing on the limit t — > oo. 

IV. SOLUTION OF THE DECOUPLED EVOLUTION EQUATIONS IN THE 

STATIONARY STATE. 

Consider the limit t — > oo and suppose that the density profiles and the stationary 
velocity of the carrier have non-trivial stationary values 

k(X) = lim k(X;t), V c = lim VJt), and A. = lim AM). (23) 

t — >oo t — >oo t — >oo 

Define next the local deviations of k(X) from the unperturbed density as 

h(X) = k(X) - Ps . (24) 

Choosing h(0) =0, we obtain the following fundamental system of equations: 

Lh(X) = 0, (25) 

which holds for A ^ {0, e±i, e±2, e±s}, while for the special sites adjacent to the carrier, i.e. 
for A = {0, e±i, e± 2 , e± 3 }, one has 

{L + A u )h{e u )+ Ps {A u -A„ u )=0, (26) 

Equations ( P5j ) and ( P6| ) determine the spatial distribution of the deviation from the unper- 
turbed density p s in the stationary state. Note also that in virtue of the symmetry relations 
in Eqs.(P) and (||), h(e ±2 ) = h{e± 3 ) and A 2 = A_ 2 = A 3 = A_ 3 . 
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The method for solving the coupled non- linear Eqs.(jj|),(f2~5|) and (^) is as follows: We 
first solve these equations supposing that the carrier stationary velocity is a given parameter, 
or, in other words, assuming that A v entering Eqs.(p5|) and ( p6|) are known. In doing so, we 
obtain h(X) in the parametrized form 

h(\) = h(\;A ±1 ,A 2 ). (27) 

Then, substituting into Eq. (|2"7|) particular values A = {e±i, e± 2 , e±3} and making use of the 
definition of in Eq.flTBD, we find a system of three linear equations with three unknowns 



of the form 

A v = l + —p u (l-p s -h(e 1/ ;A± 1 ,A 2 )), (28) 

where v = {±1,2}, which will allow us to define all A u explicitly (and hence, all h(e u )). 
Finally, substituting the results into Eq.(|5|), which can be written down in terms of A v as 

V c= ^(A-A-t), (29) 
we arrive at a closed-form equation determining implicitly the stationary velocity. 



A. Formal expression for the density profiles in the dynamic percolative environment 

as seen from the stationary moving carrier. 



The general solution of Eqs. (|25|) and can be most conveniently obtained by intro- 
ducing the generating function 

H(w 1 ,w 2 ,w 3 ) = h(\)w?w?w?, (30) 

na,n2,ri3 

where ni,n 2 and are the components of the vector A, A = e\U\ + e*2^2 + 6371,3. Multiplying 
both sides of Eqs.(^) and (|26|) by w^w^w^ 3 and performing summation, we find then that 
H(wi,W2,w 3 ) is given explicitly by 

, Y, v U v {w^-l)h{Z v )+ Ps {A v -A„ v )w^ 

H (wi, w 2 , w 3 ) = — I jr-, -, (31) 

lE,A(w^-l)-6(f + 9 ) 
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an expression which allows us to determine the stationary density profiles as seen from the 
carrier which moves with a constant velocity V c . 

Inversion of the generating function defined by Eq.(|31~D yields then, after rather lenghty 
but straightforward calculations, the following explicit result for the local deviation from 
the unperturbed density: 

h{\)=a- l {Y J A v h&)^-v- 

V 

- p s {A x - A_i)(Vi - V_i)}F(X), (32) 

where F(\) is given by 



Ai J Jo \ a 

x I„ 2 (2^aA I n3 (2^) dx, (33) 



a / \ at 



and 



a = 2_^A v -\ 



V 



= A 1 + A_ 1 +4A 2 + 6{f + 9) (34) 

Consequently, the particles density distribution as seen from the carrier moving with a 
constant velocity V c obeys 

fc(A) = p s + a~ l { J2 A u h(e u V- u - 

V 

- p s (A 1 - A_i)(Vi - V_a)}F(A), (35) 

where we have to determine three yet unknown parameters Ai, A_i and A 2 . 

To determine these parameters, we set in Eq.(|32|) A = e\, A = e_i and A = e*2, which 
results in the system of three closed-form equations determining the unknown functions A v , 
u = {±l,2}, 

6t* { , , . det C v ] , . 

A v = l + — Pv \l-Ps- Ps(A 1 - A-i)——£ , (36) 
<t det C I 



where C is a square matrix of the third order defined as 

/A 1 V„ l F(e 1 )-a A_ 1 V 1 F(e 1 ) A 2 V_ 2 F(e x ) \ 

A 1 V_ 1 F(e_ 1 ) A_iViF(e_i) - a A 2 V_ 2 F(eL 1 ) (37) 
V AiV_iF(^) A-iViF^) A 2 V- 2 F(e 2 )-a/ 

while stands for the matrix obtained from (7 by replacing the v-th column by a column 
vector ((Vi — V -\)F(e v )) v . Equation (|35|), together with the definition of the coefficients 
A u , constitutes the first general result of our analysis defining the density distribution in the 
percolative environment under study. 



B. General force-velocity relation. 



Substituting Eqs.(^) and (|37| ) into (p9|) , we find that the stationary velocity of the 
carrier particle is defined implicitly as the solution of equation: 

6r* pi det C\ - p_i det C-i 1 1 



It- 



er 



(Pi -p-i)(l-ps) U +Ps j 
t It 



detC 



where C% and C_i are the following square matrices of the third order: 

/ (Vi - V_i)F(ei) A_ 1 V 1 F(e 1 ) A 2 V- 2 F(e 1 ) \ 



(Vi - V_i)F(e_i) A_iViF(e_i) -a A 2 V_ 2 F(eLi) 
V (Vi - V_i)F(e 2 ) A_ 1 V 1 F(e 2 ) A 2 V_ 2 F(e 2 ) -a J 



and 



C_i 



A 2 V_ 2 F(e 1 ) \ 
A 2 V_ 2 F(e_!) 
A 2 V_ 2 F(e 2 ) - a/ 



(38) 



(39) 



(40) 



/ AiV_iF(ei) - a (Vx - V_i)F(ei) 
AiV-iF(gLi) (Vi-V-O^e-i) 
V A x V_iF(eb) (Vi - V_i)F(e 2 ) 
Equation ( p8|) represents our second principal result defining the force-velocity relation in 
the dynamic percolative environment for an arbitrary field and arbitrary rates of the diffusive 
and renewal processes. 
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V. CARRIER VELOCITY IN THE LIMIT OF SMALL APPLIED FIELD E, 
FRICTION COEFFICIENT AND CARRIER DIFFUSIVITY IN DYNAMIC 

PERCOLATIVE ENVIRONMENT. 

We consider now the case when the applied external field E is small. Expanding the 
transition probabilities p± and p_i in the Taylor series up to the first order in powers of the 
external field, i.e. 




(41) 



we find that V c defined by Eq. (|29f) follows 




(42) 



On the other hand, Eq.([32|) entails that 



/i(ei) - h(e-x) 



2ap s (l - p s )r* 



(3E + 0[e 2 ) 



(43) 



/r(a £(2A /a ) - A ) + 2p s 



T 



where 



A = lim E _ Q A u = 1 + — (1 - p s ), 



(44) 



and 




(45) 



while 




{P(0;3x) -P(2ei;3x)} 



(46) 



P(f; ^) being the generating function, 



P(f;0 = E^W, 



(47) 
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of the probability Pj (r) that a walker starting at the origin and performing a Polya random 
walk on the sites of a three-dimensional cubic lattice will arrive on the j-th step to the site 
with the lattice vector f || . 

Consequently, we find that in the limit of a small applied field E the force- velocity 
relation in Eq. (P8|) attains the physically meaningful form of the Stokes formula E = (V c , 
which signifies that the frictional force exerted on the carrier by the environment particles 
is viscous. The effective friction coefficient C is the sum of two terms, 

C = Co + (coop (48) 

where the first term represents a mean-field-type result Co = 6r//5cx 2 (l — p s ) (see Eq.(7)), 
while the second one, (coop, obeys 

C = l ^Ml (49 ) 

CO ° P (3a%l-p s ).(a C(2A /a )-A ) 

The second contribution has a more complicated origin and is associated with the cooperative 

behavior - formation of a inhomogeneous stationary particle distribution around the carrier 

moving with constant velocity V c . Needless to say, such an effect can not be observed within 

the framework of previous models of dynamic percolation, since there the carrier does not 

influence the host medium dynamics |l3 HT7|p!9H 22[| . 

Let us now compare the relative importance of two contributions, i.e. Co and (coop, to the 
overall friction. In Fig.2 we plot the ratio C/Co versus the creation rate / for three different 
values of the density p s , p s = 0.9, 0.7 and 0.5, while the annihilation rate is prescribed by the 
relation g = /(l — p s ) / p s . This figure shows that the cooperative behavior clearly dominates 
at small and moderate / (which entails also small values of g), while for larger /, when (/(q 
tends to 1, the mean- field behavior becomes most important. The cooperative behavior also 
appears to be more pronounced at larger densities p s . 

Consider next some analytical estimates. We start with the situation, in which diffusion 
of the environment particles is suppressed, i.e. when 1 = 0. In this case, we get 



(coop _ 2p s 



Co (1 _ p.) -l) 
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(50) 



where 



1 

y = o 




r (f + g) 
r*{l-Ps) 



) 



-i 



(51) 



Suppose first that p s is small, p s <C 1. Then, y l/3(l + r/r*(/ + p)) and we can distinguish 
between two situations: when r <C (/ + g)/r*, i.e. when the carrier moves faster than the 
environment re-organizes itself, and and the opposite limit, r ^> (/ + g)/r*, when the 
environment changes very rapidly compared to the motion of the carrier. In the former case 
we find that y w 1/3, which yields C coop /Co ~ 2p s /(6£(l/3) - 1), £(1/3) « 0.7942, while in 
the latter case we have y ~ r*/3r(/ + g>) and (coop/ Co ~ PsT*/3r(f + gr). Note, that in both 
cases the ratio (coop/ Co appears to be small, which signifies that at small densities p s the 
mean-field friction dominates. Such a result is consistent with the behavior depicted in Fig. 2 
and is not counterintuitive, of course, since in the absence of the particles' diffusion, which 
couples effectively the density evolution at different lattice sites, no significant cooperative 
behavior can emerge at small densities. On the other hand, at relatively high densities 
p s ~ 1 and t/ (1 — ps) ^> r* / (/ + g) >> r, when the carrier moves at much faster rate than 
the host medium reorganizes itself, we find that Ccoop/Co ~ T*/3r(f + g) ^> 1. This result 
stems from the circumstance that in sufficiently dense environments modeled by dynamic 
percolation a highly inhomogeneous density profile emerges even in the absence of particles 
diffusion; Here, on the one hand, the carrier perturbs significantly the particle density in its 
immediate vicinity. On the other hand, the density perturbance created by the carrier does 
not shift the global balance between creation and annihilation events, i.e. the mean particle 
density still equals p s , as we set out to show in what follows. The latter constraint induces 
then appearence of essential correlations in particles distribution and hence, appearence of 
cooperative behavior. 

Let us consider the opposite case when the renewal processes are not allowed, which 
means that the particles number is conserved and local density in the percolative environ- 
ment evolves only due to particles diffusion. In this case we find 
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Ccoop _ p s (K0\ 

Co ~~ (Zr + T*(l-p s ))(6£(l/3)-l) 



Here, the ratio Ccoop/ Co can be large and the " cooperative" friction dominates the mean- field 
one when It <C r*(3p s — 1), which happens, namely, at sufficiently high densities and in the 
limit when the carrier moves at a much faster rate than the environment reorganizes itself. 
Otherwise, the mean-field friction prevails. 

To estimate the carrier particle diffusion coefficient D c we assume the validity of the 
Einstein relation, i.e. f3D c = We find that, in the general case, the carrier diffusion 
coefficient D c reads 

Dc = ^) {, _ *jf {aoCi2Ao/ao) _ t + I^il)-} ,53, 

In the particular case of conserved particles number, when /, g — > but their ratio f/g is 
kept fixed, f/g = p s /(l — p s ), the latter equation reduces to the classical result 

D»« = {l - ^(M*C(l/S) - 1 + ^f^Y 1 } , (54) 

obtained earlier in Refs. [^] and by different analytical techniques. The result in 



Eq.(p3p is known to be exact in the limits p s C 1 and p s ~ 1, and serves as a very good 
approximation for the self-diffusion coefficient in hard-core lattice gases of arbitrary density 

It seems also interesting to analyse how random annihilation and creation of particles 
can modify the self-diffusion coefficient compared to the situation when the particles number 
is conserved. In Figure 3 we plot the ratio D^ K / D c (= C/Cnk) versus the creation rate / 
for three different values of the density p s , p s = 0.9,0.7 and 0.5. Again, the value of the 
annihilation rate g is prescribed by the relation g = f(l — p s )/p s . Figure 3 shows that the 
renewal processes affect considerably the friction coefficient and the ratio C/Cnk deviates 
strongly from the unity with the growth of the creation rate. The overall friction also falls 
off when the density increases. 

Finally, in the absence of particle diffusion (fluctuating-site percolation), our result for 
the carrier particle diffusion coefficient reduces to 
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Df r = ff2(1 6r Ps) {l - 2p.(4[(l - p.) + (/ + g)r/T*)C(y) + 3 Ps - l) -1 } (55) 

Note, however, that this result only applies when both / and g are larger than zero, such 
that the renewal processes take place. In fact, the underlying decoupling scheme is only 
plausible in this case. Similarly to the approximate theories in Refs. pl| and p5| , our 
approach predicts that in the absence of the renewal processes D^ er vanishes only when 
p s — > 1, which is an incorrect behavior. 

VI. ASYMPTOTIC BEHAVIOR OF THE DENSITY PROFILES AT LARGE 
DISTANCES IN FRONT OF AND PAST THE CARRIER. 

The density profiles at large separations in front of and past the carrier can be readily 
deduced from the asymptotical behavior of the following generating function 

+00 

N( Wl )= £ /i(n 1 ,n a = 0,n 3 = 0K 1 . (56) 

ni=— 00 

Inversion of Eq . (|31~1) with respect to the symmetric coordinates 712 and yields then 

(Axhiex) + Ps (A 1 - (wx - l) + (A- X h(€- X ) - p s (A x - A- X )) (V 1 - 1 



N(wi) 



a — A x w x — A_ x w x 

x / exp[-x]lg( 2 _f 2 x)dx + j^r^—A 

Jo a — AiW]_ — A_ x w x a — A x w x — A_iWi 

2Ao , , , 2A-. 



x 



x / exp -x l ( - — n 1 a;)fli( — - x 

Jo a — A x w x — A_ x w x v a — AiW 1 — A_ x w x 



x) — 



2A 



-I ( "J* x))dz (57) 

We notice now that N(wi) is a holomorphic function in the region Wi < u>i < W2, where 



and 
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As a consequence, the asymptotic behavior of h(nx,ri2 = 0,ri3 = 0) in the limit ni — > oo 
(resp. ni — ► — oo) is controlled by the behavior of N(wx) in the vicinity of W\ = W 2 (resp. 
Wi = Wi) (see, for example, the analysis of the generating function singularities developed 
in Ref. f3§). 



i. Asymptotics of the density profiles at large separations in front of the carrier. 



Consider first the asymptotic behavior of the density distribution of the "environment" 
particles at large separations in front of the carrier. Using the fact that 



exp 



— x 



l {yx){l 1 {x) - l (x))dx 



exp 



is a regular function when y — > 1/2, while 

I 

we find that 

(A 1 h(e 1 )+ Ps (Ai-^i))(w 2 -l 



ll(yx)dx -> - In (-— 

7T v 1 — 



2y 



N(wi) 



+ 



+ 



A^h(e-i) ~ Ps{A x - A- X ) )[m l ~ 1 



In (W2 - wi 



Then, (cf, Flajolet et al., Ref. P5[), we obtain the following asymptotical result 



/i(ni,0,0) ~ ni _ 



IT 



ni 



-ni/A-| 



where 



A + = In" 1 



a/2 - 2A 2 



+ 



V 



jCx—^w-^Tx 



'a/2-2A 2 \ Ax 



and 



+ 



47tA 2 

(A-ifrg-i) - p s (Ax - (V^ 1 - 1 

471-719 



+ 



> 0, 



(60) 



(61) 



(62) 



(63) 



(64) 



(65) 
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which signifies that the density of the "environment" particles in front of the carrier is 
higher than the average value p s and approaches p s at large separations from the carrier as 
an exponential function of the distance. 



2. Asymptotics of the density profiles at large separations past the carrier. 



We consider next the asymptotic behavior of the "environment" particles density profiles 
past the carrier particle, which turns out to be very different depending on whether the 
dynamics of the percolative environment obeys the strict conservation of the "environment" 
particles number or not (the renewal processes are suppressed or allowed). The sketch of 
this behavior is presented in Fig.4. 

a. N on- conserved particles number. In the case when partciles may disappear and re- 
appear on the lattice, one has that the root Wi < 1. We find then, following essentially the 
same lines as in the previous subsection, that 

.(A 1 h(e 1 )+p s (A 1 - A_0)(Wi - l) 



N( Wl ) 



wi-tWi 



4ttA 9 



+ 



+ 



AttAo 



111 



(66) 



Hence, in the non-conserved case the approach to the unperturbed value p s is also exponential 
when rii — > — oo, and follows 

K- 



where 



A_ = In 



7*711 ,o,o 



_! ( a /2 - 2A 2 
A., 



-|ni|/A_ 



(67) 



'a/2-2A 2 



A. 



(68) 



and 



K~ 



(A 1 h{ei)+Ps{A 1 -A^)){w 1 -l) 



+ 



+ 



4nA 2 

(A^hie-i) ~ Ps(A 1 - (Wr 1 - l) 



AuA 2 
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< 



(69) 



which implies that the particles density past the carrier is lower than the average. Note 
that, in the general case, A + < A_, which means that the depleted region past the carrier is 
more extended in space than the traffic-jam-like region in front of the carrier. The density 
profiles are therefore asymmetric with respect to the origin, n\ = 0. Since creation of parti- 
cles is favored (suppressed) in depleted (jammed) regions, while annihilation is suppressed 
(favored), one might expect that this will shift the overall density in the system, i.e. the 
average density of the "environment" particles will differ from p s . Interestingly, the overall 
deviation, i.e. the sum of local deviations over the volume of the system, of the density of 
the "environment" particles from the average value p s , appears to be equal exactly to zero, 

H(w 1 = l,w 2 = l,w 3 = l) = Q, (70) 

and hence, the driven carrier does not perturb the global balance between creation and 
annihilation of the "environment" particles. This is not, however, an a priori evident result 
in view of the asymmetry of the density profiles. 

b. Conserved particles number. Finally, we turn to the analysis of the shape of the 
density profiles of the percolative environment past the carrier in the particular limit when 
the host medium evolves only due to diffusion, while creation and annihilation of particles 
are completely suppressed. In this case, in which the particles number is explicitly conserved, 
one has that for arbitrary value of the field and particles' average density, the root Wi = 1 
and, consequently, the form of the generating function is qualitatively different from that in 
Eqs.fl6|) and 



N( Wl ) 



+ 



AttA 2 

2 Ps (A 1 -A_ 1 ' 



+ 



47rA 



Equation ( ]7T| ) implies that in the limit when the particle number is conserved the large-nx 
asymptotic behavior of /i ni) o,o is described by an algebraic function of ni with a logarithmic 
correction; that is, 
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if-ln(H) 

n>m,o,o ~ 2 ' I'^v 

where X_ is an r^-independent constant. Remarkably, the power-law decay of correlations 
implies existence of a quasi-long-range order in the percolative environment past the carrier. 
In the conserved case the mixing of the three-dimensional percolative environment is not 
very efficient and there are considerable memory effects - the host medium remembers the 
passage of the carrier on large space and time scales. 

VII. CONCLUSIONS 

To conclude, we have studied analytically the dynamics of a carrier driven by an external 
field E in a three-dimensional environment modeled by dynamic percolation on cubic lattice 
partially filled with mobile, hard-core "environment" particles which can spontaneously 
disappear and reappear (renewal processes) in the system with some prescribed rates. Our 
analytical approach has been based on the master equation, describing the time evolution of 
the system, which has allowed us to evaluate a system of coupled dynamical equations for the 
carrier velocity and a hierarchy of correlation functions. To solve these coupled equations, we 
have invoked an approximate closure scheme based on the decomposition of the third-order 
correlation functions into a product of pairwise correlations, which has been first introduced 



in Ref. for a related model of a driven carrier dynamics in a one-dimensional lattice 
gas with conserved particles number. Within the framework of this approximation, we have 
derived a system of coupled, discrete-space equations describing evolution of the density 
profiles of the environment, as seen from the moving carrier, and its velocity V c . We have 
shown that V c depends on the density of the "environment" particles in front of and past 
the carrier. Both densities depend on the magnitude of the velocity, as well as on the rate of 
the renewal and diffusive processes. As a consequence of such a non-linear coupling, in the 
general case, (i.e. for an arbitrary driving field and arbitrary rates of renewal and diffusive 
processes), V c has been found only implicitly, as the solution of a non-linear equation relating 
its value to the system parameters. This equation, which defines the force- velocity relation 
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for the dynamic percolation under study, simplifies considerably in the limit of small applied 
field E. We find that in this limit it attains the physically meaningful form of the Stokes 
formula, which implies, in particular, that the frictional force exerted on the carrier by the 
environment modeled by dynamic percolation is viscous. In this limit, the carrier velocity 
and the friction coefficient are calculated explicitly. In addition, we determine the self- 
diffusion coefficient of the carrier in the absence of the field and show that it reduces to the 
well-know result of Refs. |24|] and [^5j in the limit when the particles number is conserved. 
Further more, we have found that the density profile around the carrier becomes strongly 
inhomogeneous: the local density of the "environment" particles in front of the carrier is 
higher than the average and approaches the average value as an exponential function of the 
distance from the carrier. On the other hand, past the carrier the local density is lower than 
the average, and depending on whether the number of particles is explicitly conserved or 
not, the local density past the carrier may tend to the average value either as an exponential 
or even as an algebraic function of the distance. The latter reveals especially strong memory 
effects and strong correlations between the particle distribution in the environment and the 
carrier position. 
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Figure Captions. 

Fig.l. A generalized model of dynamic percolation. Grey spheres denote the hard-core 
"environment" particles, which perform symmetric random hopping among the sites of a 
simple cubic lattice, and can be spontaneously annihilated and created. The lighter sphere 

— * 

is the carrier, which performs a biased random walk due to an external field E, consrained 
by hard-core exclusion with the "environment" particles. 

Fig. 2. The ratio of the overall friction coefficient and the mean-field friction versus the 
creation rate for three different values of the mean density p s . The upper curve corresponds 
to p s = 0.9, the intermediate - to p s — 0.7, and the lower - to p s — 0.5. 

Fig. 3. The ratio of the overall friction coefficient and the friction coefficient in the 
conserved particles number case versus the creation rate for three different values of the 
mean density p s . The upper curve corresponds to p s = 0.5, the intermediate - to p s — 0.7, 
and the lower - to p s — 0.9. 

Fig. 4. A sketch of the asymptotic density profiles in front of and past the stationary 
moving carrier. The abbreviation CPN stands for the "conserved particles number". The 
two solid lines in the A > and A < domains denote exponential profiles, Eqs.(63) and 
(67). The dashed line in the domain A < stands for the algebraic law, Eq.(72). 
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